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ABSTRACT 

In this letter we report the detection of a strong and extremely blueshifted X- 
ray absorption feature in the 1 Ms Chandra spectrum of CXO CDFS J033260.0- 
274748, a quasar ai z = 2.579 with L2_iofceV ~ 4 x 10^^ ergs s^^. The broad 
absorption feature at ~ 6.3 keV in the observed frame can be fitted either as an 
absorption edge at 20.9 keV or as a broad absorption line at 22.2 keV rest frame. 
The absorber has to be extremely ionized with an ionization parameter ^ ~ 10^, 
and a high column density Nh > 5 x 10^^ cm~^. We reject the possibility of a 
statistical or instrumental artifact. The most likely interpretation is an extremely 
blueshifted broad absorption line or absorption edge, due to H or He-like iron in 
a relativistic jet-like outflow with bulk velocity of ~ 0.7 — 0.8c. Similar relativistic 
outflows have been reported in the X-ray spectra of several other AGNs in the 
past few years. 

Subject headings: galaxies: active — quasars: individual (CXO CDFS J033260.0- 
274748) — X-rays: galaxies 
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1. Introduction 

The 1 Ms exposure of Chandra Deep Field South (CDFS, e.g. Giacconi et al. 2002; 
Rosati et al. 2002), along with the 2 Ms exposure of Chandra Deep Field North (CDF-N, 
e.g. Brandt et al. 2001; Alexander et al. 2003) are the deepest X-ray image ever taken. With 
these two deep surveys, the cosmic X-ray background discovered by Giacconi et al. (1962), 
the origin of which has been a major goal of X-ray astronomy for almost four decades, 
is now almost completely resolved into individual sources. Most of these X-ray sources 
are extragalactic, harboring supermassive black holes. The major goal now is studying the 
properties of these X-ray sources and understanding their physical nature. 

Most recent Chandra and XMM observations have detected relativistic outflows with 
bulk velocity of Vont ~ 0.6 — 0.7c in several active galactic nuclei (Yaqoob et al. 1998, 
1999; Wang et al. 2003), revealed by extremely blueshifted emission lines due to iron or 

other elements in the X-ray spectra. In this letter, we report the discovery of a strong and 
blueshifted X-ray absorption feature (suggesting an even higher fout ~ 0.7 — 0.8c) in one of 
the GDF-S sources, GXO GDFS J033260.0-274748 (hereafter GDFSll). 

2. The Data and X— ray spectral fitting 

The 1 Ms Chandra exposure on GDF-S was composed of eleven individual AGIS obser- 
vations obtained from October 1999 to December 2000. Giacconi et al. (2002) presented the 
detailed X-ray data reduction and the final X-ray catalog (see also Alexander et al. 2003). 
The source GDFSll, ^ T from the center of the field, was observed in all the eleven expo- 
sures. Its X-ray radial intensity profile is consistent with that of a point-source (Giacconi 
et al. 2002), and no other X-ray source was detected within 45". The optical counterpart 
of GDFSll {Rvega — 21.8) was selected using the deep optical image {Rvega < 26), obtained 
with the FORSl camera on the ANTU telescope (UT-1 at VLT). Giacconi et al. (2002) 
presented the i?-band image cut-out, overplotted with X-ray flux contours (see their Fig. 
13r). We can clearly see a single point-like optical counterpart located right at the center of 
the X-ray contours, with no nearby optical source within 7". The source was firmly classified 
as a quasar ai z = 2.579 by the follow-up spectroscopy observations (see Fig. 6 of Szokoly et 
al. 2004). The X-ray to optical fiux ratio log(/2-iofeei///R) is 0.3, typical for AGNs (see, e.g. 
Hornschemeier et al. 2001). A weak, flat-spectrum radio counterpart has also been detected 
with the VLA with flux densities of 36 ±10 micro Jy at 6cm, and 33±11 at 22cm (Kellermann 
et al. in preparation), making it radio quiet with radio-to-optical ratio fiAGHz/ fn ~ 6. 

We extract the Chandra AGIS-I X-ray spectrum of GDFSll from a circle with radius of 
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7", which is the 95% encircled-energy radius of the ACIS point-spread function at the source 
position. The local background was extracted from an annulus with outer radius of 19" and 
inner radius of 9". In Fig. 1 we present the summed spectrum (source plus background) and 
the background evaluated in the outer annulus. The source is fairly bright in X-ray with 
~ 1040 net X-ray counts in the 0.5 — 2.0 keV band and ~ 350 in the 2.0 — 9.0 kcV band, 
allowing us to perform X-ray spectroscopy. During the fit, we use the C statistics (Cash 
1979; Nousek & Shue 1989), which has better performance with respect to the analysis, 
particularly for spectra with low detected counts. We generate the X-ray telescope response 
and ACIS-I instrument response for each single Chandra observation, and sum the response 
files weighting them for the corresponding exposure times. The final time-weighted response 
files were used for spectral analysis. We utilize XSPEC vll.2 to perform the spectral fitting. 
All the spectral fitting was done in the energy band 0.5 — 9.0 keV and all the statistical 
errors quoted in this paper are at the 90% confidence level for one interesting parameter. 

We first fitted the spectrum with a simple power-law plus a neutral absorber in the 
source frame. A Galactic neutral hydrogen absorption column of 8 x 10^^ cm~^ (Dickey & 
Lockman 1990) was also included. The results are presented in Table 1. The spectrum was 
well fitted by a power-law (T = I.TIq}) with weak absorption {Nfj = O.I^'^q'i' x 10^^ cm~^). 
The intrinsic, rest-frame 2.0 — 10.0 keV luminosity is 4.3 x 10^^ ergs s^^ {Hq = 70.0 km 
s~-'^Mpc~^, Q.rn = 0.27, Q.A = 0.73). The best-fit continuum model and the ratio of data to 
model is shown in Fig. 2. Both in Fig 1 and 2, we can see a strong absorption feature at ~ 
6.3 keV in the observed frame. The absorption appears to be optically thick, since the ratio 
of data to model reaches zero at energies around 6.3 keV. 

We add an absorption edge to our fit, which attenuates the continuum above Eedge 
with the optical depth r = To{E / Eedge)~^ ■ An absorption edge with Eg^ge — 5.8 keV and 
To = 3.5 significantly improves the fit (AC = — 15 with two extra free parameters, see Table 
1). We also tried to model the absorption feature by a saturated absorption line model. The 
model notch of XSPEC was used by fixing the covering fraction at 0.99 to represent a blank 
absorption trough. The fitting is slightly better than the edge model with AC = — 17 for 
two extra free parameters, which are located at E^ = 6.3 keV, and the line width of 0.9 keV 
(both in the observed frame). The line width can also be taken as the equivalent width of 
the absorption line since the absorption is saturated. We also perform spectral fits to search 
for possible other absorption edges/lines at lower energies. We try different energy entries 
from 0.5 to 6 keV, but we find that no further absorption edge/line is statistically required 
(AC < 3). 

The X-ray spectrum of CDFSll was found to be variable with high probability (> 3(7, 
Paolillo et al. 2004). However, due to the limited number of X-ray photons over 5 keV, 
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and the fact that the absorption feature is optically thick with black trough, we are unable 
to study the variability of the absorption itself, based on our data. Sometimes, a broad 
emission line might actually be mimicked by a strong absorption edge at higher energy (e.g.. 
Reeves, Porquet & Turner 2004) and vice versa, especially in X-ray spectra with low S/N 
or low energy resolution. To check if this is the case for CDFSll, we fit the spectrum with a 
power law continuum and a broad emission line. We tried different central values for the line 
energy from 4 to 6 keV, and found that an emission line can only improve the C statistics 
by AC < —2. This indicates that the absorption feature in CDFSll can't be ascribed to 
the presence of a broad emission line. 



3. Discussion 

We discuss here the possible origin of the significant broad absorption feature we de- 
tected in the X-ray spectrum of CDFSll with the IMs Chandra ACTS exposure. The 
feature locates at energies > 20 keV in the rest frame, with a confidence level of 99.98% 
according to F-test. We note that the confidence level given by F-test for an absorption 
feature might not be accurate (Protassov et al. 2002). Here we re-estimate the confidence 
level of the absorption feature via a simple approach. In the energy range 5.8 — 6.8 kcV 
where the absorption feature is located, we detected a total of 11 photons (source + back- 
ground), while the continuum model plus background predict 32 photons. The cumulative 
Poisson probability for such a deviation is ~ 1.7 x 10~^. Since the absorption feature was 
examined over the whole spectral band, which is 0.5 — 9.0 keV, the probability of detect- 
ing such a broad spurious absorption feature randomly in the spectrum is approximately 
1.7 X 10~^/(9.0 - 0.5) = 2.0 X 10~^. This is consistent with the confidence level obtained 
with the F-test. 

We also perform extensive Monte Carlo simulations to check the confidence level. First, 
we simulate 10^ artificial spectra based on the best-fitting continuum model. We then 
search for spurious absorption lines in the artificial spectra by adding a notch component. 
All the three parameters of the notch model (central energy, width, and covering factor) are 
thawed. We did not fix the covering factor at 0.99 during the simulation because shallower 
absorption line (i.e., with lower covering factor) can reach comparable significance level at 
lower energy due to the higher S/N. To ensure an efficient search over the whole band, we 
perform the search in narrow energy bins (such as 0.5 — 1.5 kcV, 1.5 — 2.5 kcV, etc) with width 
comparable to that of the observed absorption feature, and count the total number of spectra 
in which we detected spurious absorption line with AC < —17. We found only 6 spectra 
with statistically significant broad (energy width > 0.5 keV) spurious absorption lines. This 
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indicates a confidence level of 99.94% for the broad absorption feature we detected. Note also 
that, while fitting the real spectra, we use only 2 free parameters of the notch model, while in 
the simulation, we use 3. We remark that the confidence level would be higher (99.98%) if wc 
searched for spurious features with AC < —19.6 for 3 free parameters, instead of AC < —17 
for 2. 

We note that there is also weak evidence of the absorption feature in the 370 ks XMM 
spectrum (Streblyanska ct al. 2004) of CDFSll. The XMM spectrum is significantly steeper 
(r = 2.3) and it was obtained about one year later. The statistical quality of the XMM 
spectrum around 6 kcV is much lower due to the high background noise and the steeper 
continuum. By fitting the Chandra and XMM spectra simultaneously, we found an even 
higher confidence level of 99.997% for the absorption feature, based on F-test. However, due 
to the very limited statistical quality, we are unable to use the XMM spectrum independently 
to constrain the nature and possible time variability of the absorption feature. 

We conclude that the confidence level of the absorption fine we detected in CDFSll is 
> 99.98%. Even considering the number of all the CDFS spectra we examined for interesting 
features (~ 30 CDFS spectra have X-ray photons > 500, whose quality is high enough to 
enable the search for broad line features), the absorption feature is still significant with a 
confidence level of 99.4%. 



3.1. Instrumental Artifact? 

Could such an unusual X-ray feature be due to any instrumental artifact or to some 
aspect of our analysis? We first consider improper background subtraction as a possible origin 
of the absorption feature. We check that we obtain consistent results when using background 

spectra extracted from different regions of the detector. Wc also note that the absorption 
feature is already significant in the spectrum without background subtraction. Finally, we 
do not see any bump at ~ 6 keV in the background spectrum cither, suggesting that the 
significance level of the feature is not magnified by background subtraction. Therefore, 
we conclude that the absorption feature is not related to uncertainties in the background 
subtraction. 

Could the absorption feature be an artifact of calibration uncertainty in the ACIS in- 
strumental response function? As far as we know, no such an artifact (i.e., an absorption 
feature at ~ 6.3 keV) has been reported. We examined the spectra of nearby X-ray sources, 
but we found no evidence of absorption feature at ~ 6.3 keV among them. Tozzi et al. (2005) 
presented the X-ray spectra of CDFS sources using updated Chandra calibration files (with 
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CIAO3.0.1 and CALDB2.26 instead of CIAO2.0.1 and CALDB2.0 adopted in Giacconi et al. 
2002 and in this letter). We repeated the fitting using the updated spectrum and obtained 
consistent results for both the continuum and the absorption feature. This confirms that the 
absorption feature is not due to improper data calibration. Note that Wang et al. (2003) 
reported a puzzling strong emission line at ~ 6.3 keV in the observed frame in CXO CDFS 
J033225.3-274219 (CDFS46). The two unusual features (the emission feature in CDFS46, 
and the absorption feature in CDFSll reported in this letter) located at similar energies in 
the observing frame, could suggest a common origin related to instrumental effects. However, 
there are strong evidences against this occurrence: since the total exposure was composed of 
eleven individual observations with different roll angles, photons from each source fall into 
different positions on different chips during different observations. It is unlikely that an un- 
known instrumental artifact (if there is any) would affect photons from all eleven exposures. 
Furthermore, the two sources are around 9' apart, and fall into different chips during the 
observations. We conclude that the similar energy of the two features is just a coincidence. 



3.2. The nature of the absorber 

Among the heavy elements, iron is the only one that can produce such a strong absorp- 
tion feature at high energies. However, the rest-frame energy of the absorption (> 20 keV) 
is too high for static iron in the rest frame, since the highest energy transition is the K-shell 
edge of H-like ion at 9.28 keV. Therefor we consider two possibilities: either the absorber 
locates at much lower redshift along the line of sight, or in an relativistic outflow intrinsic 
to the quasar. 

The continuum fitting yields a marginal intrinsic absorption with Nh ~ 10^^ cm~^, 
which is too small to account for the strong absorption feature at ~ 6.3 keV. We can't 
locate other absorption features either in the spectrum at lower energy (due to relatively 
lower ionized ions, such as the OVHI absorption edge at 0.87 keV). This suggests that the 
absorber has to be extremely ionized, with almost all abundant elements fully ionized, and 
the absorption is dominated by H or He-like iron atoms. Adopting the photon ionization 
absorption model absori (Magdziarz & Zdziarski 1995) to fit the spectrum, we obtain an 
ionization parameter^ ^ > 9000 and Nh > 3 x 10^^ cm~^ (Table 1). Note the upper limits 
of ^ and Nh are poorly constrained since the two parameters are degenerate. Assuming 
^ = 10000, we obtain Nh = 5 x 10^"^ cm~^. We assume a temperature of 10^ K and 



= L/nE?, where L is the integrated incident luminosity between 5 eV and 300 keV, n is the density 
of the material and R is the distance of the material from the illuminating source (Done et al. 1992). 
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solar abundance for the absorber. We note that an intrinsic absorption at the level of 
Nh = 5 X 10^^ cm~^ (well within the Compton-thick regime) would strongly attenuate the 
continuum photons. However, if the iron in the absorber is a factor of 10 overabundant with 
respect to the solar value, we obtain Nh ~ 5 x 10^^ cm~^. 

The absorption feature can also be fitted by a broad saturated absorption line. The 
most likely responsible lines are 6.70 Fcxxy Ka (6.70 keV) and/or Fexxy/ Kci (6.97 keV). 
Assuming an average depth of the absorption line r = 2, and a factor of 10 overabundant for 
iron, we obtain Nh ~ 5 x 10^'^ cm^^. We point out that for such an absorber, we also expect 
Fcxxy K/3 (7.85 keV), 'Pexxvi K/? (8.26 kcV) absorption lines, and strong Ycxxv (8.83 
keV), Yexxvi (9.28 keV) absorption edges. Given the few counts we measure above 7 keV, 
we can't check the existence of these features. Future X-ray missions with higher sensitivity 
can help us to unveil the nature of the absorber by showing us the spectral properties at 
energies > 7 keV and the absorption profile with higher S/N. 

If the absorption is due to foreground absorber, the absorber has to be located at a 
much lower redshift ~ 0.5 (if due to H or He-like iron absorption edges) or ~ 0.1 (if due to 
H or He-like iron resonant absorption line). The absorber must be extremely ionized with 
Nh > 5 X 10^^ cm"^. Note the foreground absorber is unlikely photo-ionized, otherwise we 
should have seen the extra photon ionization source, thus it has to be extremely hot with 
a temperature of ~ 10^ K to reach the required ionization stage. Such a high temperature 
and Nh is very unusual for intervening systems. The only possible candidate is the hot 
intra-clustcr medium found in clusters of galaxies (sec Rosati, Borgani & Norman 2002 and 
references therein). However, for a cluster with T = 10^ K, we expect an X-ray luminosity of 
~ 3 X 10^^ erg s"^ (based on the Lx — T relation, e.g. Lumb et al. 2004). If located at z ~ 0.5, 
such a cluster would be around 300 times brighter than CDFSll, and its emission would be 
obviously prominent in the Chandra image. We conclude that the X-ray absorption feature 
we detected in CDFSll can't be due to an intervening systems along the line of sight. 

Finally, we consider a relativistic outflow intrinsic to the quasar, with speed of ~ 0.7c 
(Doppler Factor DF — 2.3 for absorption edge) or ~ 0.83c {DF — 3.2 for absorption line), 
to produce the observed blueshift of the absorption feature. Note that similar blueshifted 

featTircs have been reported in the X-ray spectra of 3 AGNs: a highly blueshifted Ovn 
emission line in PKS 0637-75 with DF ~ 2.7 - 2.8 (Yaqoob ct al. 1998); a blueshifted F(^K 
emission line {DF ~ 2.4 - 2.6) in QSO PKS 2149-306 (Yaqoob et al. 1999); and a strong 
blueshifted Fc-K emission line {DF ~ 2.3 - 2.5) in GXO GDFS J033225.3-274219 (Wang 
et al. 2003). The good agreements of the blueshift factors from different sources strongly 
suggest a close origin of the relativistic outflow, and strengthen the statistical significance of 
these features. 
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Taking log^ = 4 and Nh = 5 x 10^^ cm~^, we estimate the distance of the absorber 
from the central source R by assuming the thickness of the absorber AR ~ O.li?. We 
obtain i? ~ 6 x 10^^ cm. Following Reeves, O'Brien & Ward (2003), we estimate the rate 
of relativistic outflow. The outflow rate (Mout) can be calculated as Afout = ^'^out^p-^x/'C- 
Assuming that the outflow subtends a solid angle of 0.1 steradian, wc obtain a outflow rate 
of ~ 10^^ g s^^ or ~ IOMq year~^, similar to that of the outflow in PDS 456 (wout ~ 0.17c, 
Reeves et al. 2003). However, because of the extremely high outflow speed in this case, 
the kinetic energy carried by the relativistic material outflow would be ~ 2OM0 year~^, and 
then requires an accretion rate of at least 200 Mq year~^ with an assumed energy production 
efficiency e = 0.1 which entirely goes into the outflow. Such an accretion rate is far beyond 
the Eddington limit even for a supermassive black hole of KfMQ (~ 3Mq ycar~^). To keep 
the rate of kinetic energy outflow lower than that released by Eddington accretion, a more 
reasonable outflow rate would be at least 100 times smaller (< O.IMq year"-*^), and so does 
the solid angle (< 0.001 steradian). The extremely high velocity and small covering factor 
strongly link the physics of the outflow to relativistic jets in AGN and X-ray binary systems 
(such as the famous XRB SS 433, see Mighari, Fender & Mendez 2002). 

We conclude that the X-ray absorption feature in CDFSll is due to an ionized jet-like 
outflow intrinsic to the quasar, with a bulk velocity of ~ 0.7 - 0.8 c. Similar outflows have 
been reported in the X-ray spectra of several other AGNs. Future X-ray mission with higher 
intensity can help us understand the nature of the outflows by providing X-ray spectra with 
higher S/N. 

The work of JW was supported by Chinese NSF through NSF 10473009 and the CAS 
"Bai Ren" project at University of Science and Technology of China. 
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Table 1. Spectral fits to CDFSll 
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^The covering factor of the notch model was fixed 
to 0.99 to represent a heavily saturated absorption. 

^The upper limit of the absorption depth was 
poorly constrained. 
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Fig. 1. — The summed (source plus background) X-ray spectrum of CDFSll and the ex- 
pected background (points with dotted error bars). The spectra were rebinned for display 
purpose. 
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Fig. 2. — The spectral data (rebinned for display purpose), best-fit continuum models and 
the ratios of data to model for CDFSll. 



